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Abstract 

In this paper we study the scattering theory of the classical hyperbolic Sutherland model as- 
sociated with the C n root system. We prove that for any values of the coupling constants the 
scattering map has a factorized form. As a byproduct of our analysis, we propose a Lax matrix for 
the rational C n Ruijsenaars-Schneider-van Diejen model with two independent coupling constants, 
thereby setting the stage to establish the duality between the hyperbolic C n Sutherland and the 
rational C n Ruijsenaars-Schneider-van Diejen models. 
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1 Introduction 



In the study of interacting many-particle systems it is hard to overestimate the importance of scattering 
theory. At the same time, it is notoriously difficult to obtain rigorous results in this subject. It is 
a very fortunate situation that for certain integrable systems, defined on the real line, the scattering 
theory has been completely understood. In particular, we have full control over the scattering theory 
of the Toda systems, the Calogero-Moser-Sutherland models, and the Ruijsenaars-Schneider models 
associated with the A n root system (see e.g. pQ, [2], [3], [JJ, [5]). Their characteristic feature is that 
the scattering map has a factorized form. Though these integrable many-particle systems have natural 
generalizations to other root systems as well (see e.g. [6]), the scattering theory of the non-^4 n -type 
models is far less developed than that of the A n -type systems. 

In this paper we undertake the task to understand the scattering behavior of the classical hyperbolic 
C n Sutherland model. (For background information on the C n -type Sutherland systems see e.g. [6], 
[Z]) [§]) [S]-) Recall that the phase space of this model is the cotangent bundle of the Weyl chamber 
c := {q = (qi, . . . , q n ) € W 1 \ q± > . . . > q n > 0} and the dynamics is governed by the Hamiltonian 

c^l i<t<b<n ^ smh ~ qb > smh ^ a + 1b)J ^ smh (2q c ) 

where g and gi are arbitrary non-zero real numbers, the so-called coupling parameters. By the 
repulsive nature of the interaction we expect that the particles move asymptotically freely for very 
large positive and negative values of time t, thus it makes sense to study the scattering map that relates 
the asymptotic phases and momenta of the past and the future. Using only elementary algebraic 
techniques, the main goal of the paper is to show that the scattering map of the C n -type model 
also has a factorized form, i.e., the classical phase shifts are entirely determined by the two-particle 
processes and by the one-particle scatterings on the external field. The precise statement is given in 
Theorem 3. Though this result does meet our expectations, to our knowledge, its rigorous proof has 
not appeared in the literature before. 

To understand the scattering properties of the hyperbolic C n Sutherland model we closely follow 
Ruijsenaars' seminal work [3] on the A n system. One of the upshots of his approach is that it reveals a 
natural action-angle duality between the hyperbolic Sutherland and the rational Ruijsenaars-Schneider 
models. Surprisingly, for root systems other than A n , relatively few is known about the duality 
between the Sutherland and the Ruijsenaars-Schneider-van Diejen (RSvD) models. Even the Lax 
representation of the generic RSvD dynamics is missing, except for some very special one-parameter 
family of C n and BC n models obtained by the natural ^-folding of the original Am—\ and Am systems 
|10| . Only partial results [11] are known for the D n root system, too. However, as a byproduct of our 
scattering theoretic analysis, we obtain a natural candidate for the Lax matrix of the two-parameter 
family of rational C n RSvD models. The Lax matrix presented in Lemma 2 generalizes the known 
one-parameter family of rational Lax matrices obtained by the folding procedure and it generates the 
rational C n RSvD Hamiltonian [12] with two independent coupling parameters. 

The paper is organized as follows. To keep the presentation self-contained, in Section 2 we collect 
the necessary background material and fix the notational conventions. In Section 3 we present our 
results on the scattering theory of the hyperbolic C n Sutherland model. After the discussion part in 
Section 4, we finish the paper with an Appendix on some useful facts from linear algebra. 
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2 Preliminaries 



In this section we gather some basic facts about the hyperbolic C n Sutherland model. We start with 
a short review on some group theoretic material related to the non-compact Lie group U(n,n), then 
we discuss the Lax representation of the Sutherland dynamics. For the details the reader may consult 

m m, ®, m, s, m- 

2.1 Group theoretic background 

Take an arbitrary n G N = {1, 2, . . .} and let N := 2n. With the aid of the unitary matrix 

C: =U o) €t/(JV) (2) 

we define the non-compact real reductive matrix Lie group 

U(n,n):={yeGL(N,C)\y*Cy = C}. (3) 

The fixed-point set of the Cartan involution 0(y) := (y -1 )* is the maximal compact subgroup 

U(n, n)+ := {U G U(n, n) \ U is unitary } S U(n) x U(n), (4) 

meanwhile the submanifold U (n, n)_ := {y G U (n, n) \ Q(y) = y^ 1 } consists of the Hermitian elements 
of U(n,n). Notice that C ([2]) is central inside U(n,n)+. 

The Lie algebra of U(n, n) has the form u(n, n) = {X £ gl(N, C) | X*C + CX = 0}. The natural 
trace-pairing (X,Y) := tr(XY) provides an invariant, non-degenerate, bilinear form on u(n,n). The 
Lie algebra involution 6{X) = —X* corresponding to induces the orthogonal ^-gradation 

u(n, n) = u(n, n) + © u(n, n)_ (5) 

with the eigenspaces u(n,n)± := ker(0=pld). The bilinear form ( , } is negative definite on the sub- 
algebra u(n, n)+ and positive definite on the complementary subspace u(n, n) Note that u(n, n) + 

(resp. u(n,n)_) consists of the anti-Hermitian (resp. Hermitian) elements of u(n, n). 
Now with any real n-tuple q = (q\, . . . , q n ) G M. n we associate the diagonal matrices 

q := diag(gi, . . . ,q n ) G s((n,R) and Q := diag(q, —q) G Ql{N, E). (6) 

The subset o := {Q = diag(g, — q) \ q G W 1 } is a maximal Abelian subspace in u(n, Its centralizer 
inside U(n, n) + is the Abelian group 

M := {diag(e ix , e ix ) | x S M n } < U(n, n)+ (7) 

with Lie algebra m := {diag(i%, i%) | x G R n } < u(n, n) + . Notice that both o and m are realized by 
diagonal matrices. Let (resp. m^) denote the subspace of the off-diagonal elements of u(n, n)_ 
(resp. u(n, n) + ), then we can write u(n, n)_ = a © a 1 - and u(n, n) + = m © m^. 

The subspace m 1 - © a 1 - formed by the off-diagonal elements of u(n, n) is invariant under the linear 
operator adQ, for any Q G o. Therefore the restricted operator 

ad Q := ad Q | m ± ffia ± (8) 
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is well-defined, with spectrum 

£7(adg) = {q a - q b , ±(q a + q b ), ±2q c \ a, 6, c G N n , a ^ b}, (9) 

where N n := {1, . . . , n} C N. The regular part of o is defined by the subset 

dreg := {Q G a | adQ is invertible } C a. (10) 

Since a \ a reg is a union of finitely many hyperplanes, a reg is open and dense in o. The subset 

c:={Q = diag(q, -q)\q= (<?i, . . . , q n ) G K n , q x > . . . > q n > 0} (11) 

is a connected component of a reg , i.e., it is an open Weyl chamber. Note that the configuration space 
of the hyperbolic C n Sutherland model (pQ) can be identified with c. In the following we will frequently 
use the identification c = {q = (q±, . . . , q n ) £W n \qi>...>q n >0}. 

As is known, the elements of u(n, n)_ can be 'diagonalized' by conjugation with elements from 
U(n,n) + . More precisely, let c~ denote the closure of c, then the map 

c" x U(n, n)+ 3 (Q, U) ^ UQU- 1 G u(n, n)_ (12) 

is well-defined and onto. Moreover, the regular part of u(n,n)_, 

(u(n,n)_) rcg := {UQU^ 1 \ Q G c, U G U(n,n) + } C u(n,n)_, (13) 

is an open and dense subset inside u(n, n)_, admitting the smooth bijective parametrization 

c x (U(n, n) + /M) 3 (Q, UM) h-> UQU' 1 G (u(n, n)_) reg . (14) 

That is, the above diffeomorphism provides the identification (u(n, n)_) reg = c x (U(n,n) + /M). 

2.2 Lax representation of the Sutherland dynamics 

Let E G C N denote the column vector with components E a = 1, E n+a = — 1 (a G N n ), and set 

i := ig(EE* - 1 N ) + i(g - g 2 )C G m ± C u(n, n) + , (15) 

where 5 and 52 are arbitrary non-zero real parameters. Utilizing the Riesz-Dunford functional calculus, 
to any (q,p) G T*c = c x M n cK"xR" we associate the N x N matrix 

£(q,p) '■= P — coth(adQ)£ Gu(n,n)_. (16) 
Consider the m- valued function <P := idiag(</?i, . . . , V9 n , </?i, . . . , <p n ) with 

(Pc{q,p) = -9 ^2 (sinh(g c -g a )~ 2 + sinh(g c + g a )" 2 ) -5 2 sinh(2g c ) -2 , (17) 

oeN„\{c} 

and define the N x N matrix 

B(q,p) := <P(q,p) + sinh(adg)~ 2 ^ G u(n, n) + . (18) 
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As is known, the matrix-valued functions C and B provide a Lax pair for the hyperbolic C n Sutherland 
model. More precisely, along a smooth regular curve q(t) 6c (t G R) the Lax equation C = [C,B] is 
satisfied if and only if q(t) is a solution of the hyperbolic C n Sutherland dynamics. 

An obvious consequence of the Lax representation of the dynamics is that the solution curves can 
be realized as projections of certain geodesies on the Riemannian manifold f/(n,n)_. Indeed, take an 
arbitrary solution q(t) € c of the Sutherland dynamics and set p{t) := q(t). The differential equation 
n(t) _1 n(t) = B(q(t),p(t)) has a unique smooth solution u(t) G U (n, n) + (t G R) with initial condition, 
say, u(0) = In- Now for any t G R we define the positive definite matrix 

y(t) := u^e^uit)' 1 G U(n,n)-. (19) 

Clearly y(t) is a smooth function of t satisfying the equation y~ l ij + yy~ l = \uLu~ x . It follows that 
y(t) is a solution of the geodesic equation on U(n,n)-, i.e., 

d fy-H + yy-^ u{ t-[c,B\)u- l = ^ (20) 



dt V 4 

with y(0) = e 2Q (°) and y(0)" 1 y(0) + y(0)y(0)" 1 = 4£(g(0),p(0)). However, upon introducing 

L(g,p) := cosh(ad Q ) _1 /:(g,p) = P - sinh(adQ) -1 ^ G u(n, n)_, (21) 
the unique solution of the geodesic equation with the above initial conditions is the curve 

y (t) = e Q(°) e 2*£(9(o), P (0)) e Q(o)_ ^22) 

Comparing ()19p and f)22[) we see that Q(t), and so the trajectory q(t), can be recovered by diagonalizing 
the matrix flow (|22p . In particular, we have the spectral identification 

{e 2qi{t \ . . . ,e 2qn{t \e~ 2qn{t \ . . . ,e~ 2qi(t) } = cr(e 2Q(0) e 2iL(<?(0) ' p(0)) ), (23) 

whence the temporal asymptotics of the trajectory q(t) can be understood by analyzing the temporal 
asymptotics of the eigenvalues of the matrix flow (|22p . Though it is a natural matrix analytic question, 
to our knowledge the first reference containing the solution of this problem is Ruijsenaars' paper [1]. 



3 Temporal asymptotics 

In this section we work out the temporal asymptotics of the hyperbolic C n Sutherland dynamics. 
Our main guide is Ruijsenaars' result on the temporal asymptotics of the eigenvalues of exponential 
matrix flows (|22p . As dictated by Theorem A2 in [3J, the plan is to find the matrix entries of e 2 ®^ 
in an orthonormal basis, in which L(q(0),p(0)) is diagonal with decreasing diagonal entries. Having 
control over the matrix entries of e 2 ®^ in this new basis, simply by computing the quotients of the 
consecutive leading principal minors, we can determine the temporal asymptotics of the eigenvalues 
of the matrix flow. 

In the rest of the paper we simply write L and Q in place of L(q(0),p(0)) and Q(0). Since the 
initial conditions q(0) and p(0) = q(0) can be arbitrary, we think of L and Q as matrix-valued smooth 
functions over the phase space T*c = c x R n . 



5 



3.1 Analyzing the spectrum of L 

Recall that it is a crucial assumption of Theorem A2 in [3] that the spectrum of matrix L is simple. 
To examine the spectrum of L, we set up an equation for L G u(n, n)_ and A := e 2< ^ £ C/(n,n)_ as 
follows. By applying the linear operator sinh(adg) on L (|2ip . we get 

sinh(ad Q )L = (24) 

which entails Le 2 ® — e 2 ®L = 2e^£e^. Now, recalling (fl~5|) . we can write 

2i 5 vl + LA — AL = 2ig(e Q E)(e Q E)* + 2i(g - g 2 )C. (25) 

Note that this equation is the complete analogue of Ruijsenaars' commutation relation (equation 
(2.4) in [3]) he analyzed to discover the remarkable action-angle duality between the hyperbolic A n 
Sutherland and the rational A n Ruijsenaars-Schneider models. If g 2 = g, the right hand side of (|25p is 
a matrix of rank one, therefore the analysis of the equation is relatively straightforward. Notice that 
the special case g 2 = g corresponds to the Z2-folding of the A 2n _\ model. However, when g 2 ^ g, the 
innocent looking term 2\{g — g 2 )C complicates the analysis considerably. 

To proceed further, we diagonalize L G u(n, n)_. As we saw in (|12p . we can write 

L = ULl/- 1 = ULU* (26) 

with some Let" and U G U{n,n) + . Observe that L is unique, having the form L = diag(A, —A) 
with some Ai > . . . > A n > 0, but the choice of U is not unique. However, at this point all we need is 
the existence of the pair (L, U), the consequences of the non-uniqueness will be discussed at the end of 
this subsection. Now from ([25]) we conclude that L and A := XJ~~ X AXJ € f/(n,n)_ satisfy the equation 

2\gA + LA - AL = 2ig{U* e Q E){U* e Q E)* + 2i{g - g 2 )C. (27) 

Computationwise it is very fortunate that the matrix A~ l E U(n, n)_ obeys a similar equation. Indeed, 
by conjugating the above equation with C, we get 

2igA- 1 - LA- 1 + A- X L = 2ig(CU*e Q E)(CU*e Q E)* + 2i(g - g 2 )C. (28) 

Upon introducing the purely imaginary numbers 

x c = —x n + c ■= (2ig)~ 1 \ c G iR (c <E N n ), (29) 

the column vector 

F := U*e Q E G (30) 

and the real parameter 

e := 1 - g^ 1 G M, (31) 
for the matrix entries of A and yl -1 we obtain 



- F k F l + eC k ,i r x {CF) k {CF) l + eC Kl 

Ak,l = — , [A ) k i = ■ , (32) 

1 + X k ~ Xl 1 - Xk + Xl 



for any k, I G Nat. Note that the relations 

JV 



F k F j + eC ktj (CF) J {CF\ + eC hl _ ^ 

Z — / 1 + X k — Xj 1 — Xj + Xl 



3=1 
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obviously follow from (|32p . 

Having equipped with the above formulae, we are able to analyze the spectrum of L and the 
properties of column vector F. For convenience, we introduce the notations 

f c := F c , h c := F n+C , z c := f c h c (c G N n ). (34) 

Lemma 1. If the non-zero coupling parameters g and g<i satisfy g<i ^ 2g, then the components of the 
column vector F are non-zero and L is a regular element o/u(n,n)_. 

Proof. We only show that the components of F are non-zero. Proving by contraposition, suppose 
that z c = for some c G N n . With k = I = c, from ([33|) we get the quadratic relation e 2 = (1 + 2x c ) 2 , 
i.e., e = ±(1 + 2x c ). By comparing the real parts we obtain e = ±1, which contradicts our assumption 
on the coupling parameters. Along the same line, by appropriately specializing the indices in equation 
(f33j) . the regularity of L also follows. ■ 
Remark. Henceforth we assume that the parameters g and g<i satisfy the additional technical condi- 
tion g2 7^ 2g. Notice, however, that it does not restrict the values of the physically relevant positive 
coupling constants g 2 and g 2 , since the pairs (g,g2) and (g,— 52) generate the same couplings in the 
model (P). In principle, without loss of generality, we could have imposed the condition gg2 < at 
the outset, thereby automatically excluding the case g2 = 2g. 

To conclude this subsection we wish to point out that the construction of z c (|34p results in a 
well-defined smooth function on the phase space T*c. As we saw in (|14p . by the regularity of L, the 
non-uniqueness of the diagonalizing matrix U G U(n, n)+ defined in ([26]) is controlled entirely by the 
centralizer subgroup M ([7]). Namely, the only freedom in the choice of U can be characterized by the 
transformations 

Un> Udmg(e ix ,e ix ), (35) 

generated by some x £ Now observe that the components f c and h c f)34|) of the column vector F 
(j30|) transform as 

f c ^ e -ixc fc and hc ^e- iXc h c , (36) 

hence z c = f c h c is independent of the choice of the representative U. It means that to each L = L(q,p) 
we can associate the non-zero complex numbers z c = z c (q,p) (c G N„) in a unique and well-defined 
manner. To show that their dependence on the phase space variables (q,p) is smooth, we notice that 
by choosing appropriate smooth local sections of the (smooth) fiber bundle 

c x U(n,n) + -» c x (U(n,n) + /M) ^ (u(n, n)_) reg , (37) 

we can work with representatives U G U(n, n) + depending smoothly on the phase space variables in 
a small neighborhood of any given (q,p) G T*c. Thus, z c is smooth around any (q,p), proving its 
smoothness over the whole phase space. 

3.2 The structure of A 

In this subsection we proceed with a detailed analysis on the structure of matrix A. To this end, 
we make use of Jacobi's theorem in linear algebra, i.e., we exploit some non-trivial relations between 
certain minors of matrices A and B := (j4 _1 ) t . For convenience, in Appendix A we provide a brief 
account on the relevant theorems from linear algebra. 
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In the following we keep the index c 6 N n fixed and apply Jacobi's theorem on appropriate minors 
of A and B. Since det(A) = det(e 2< ^) = 1, let us keep in mind that det(^4) = 1 also holds. To make 
the presentation shorter we introduce the notations 

n n n , , 

■H \\ l + x a -x b {l + x c -x a )(l+x c + x a ) 

d=l a, 6=1 a=l 

Clearly T> c € R \ {0} and w c € C \ {0}. Note also that in the equations below the matrices ej-,i stand 
for the elementary matrices, i.e., their entries are defined as (ekl)k' V = $k,k'fii,i'- 
As a first application of Jacobi's theorem, we can write 

b( j l n )=-A( n t) U + C l H ). (39) 

V 1 ••• n + c ••• n J y n + 1 ••• c ••• 2n J 

Let R and S denote the n x n submatrices corresponding to the above minors on the left and on the 
right, respectively. Upon introducing the Cauchy-type n x n matrix ^ with entries 

* a , b := : ^ hh _ ifMc, and * , c := — , (40) 

from (|32p we see that 

^ = ^ + e(l + 2x c )- 1 e CjC , (41) 

meanwhile the entries of S can be identified as S a ^ = R a ^ (a, b € N n ). Therefore, equation (|3"9"j) can 
be cast into the particularly simple form 



det(-R) + det(i?) = 0. (42) 

Since R is a rank one perturbation of the determinant formula (|A.6p spells out as 

det(fl) = det(^) + e(l + 2x c )~ 1 C CjC , (43) 

where C CiC is the cofactor of ^ associated with entry ^ C)C . Since \& is of Cauchy-type, we get 

det(*) = (1 + 2x c )~ 1 V c oj c z c and C c . c = £> c . (44) 

Plugging these formulae into (|42|) . for z c we obtain the linear relation 

(1 - 2x c )uj c z c + (1 + 2x c )lJ c z c + 2e = 0. (45) 

Notice that this single equation does not determine uniquely the complex quantity z c . 

To get an independent relation for z c , we turn to Jacobi's theorem, again. Namely, we can write 

b{\ n n + c )=A( n + 1 55 2n V (46) 

y 1 ••• n n + c J \ n + 1 ••• n + c ■■■ 2n J 

where the symbol n + c means that the indicated row (and column) is deleted in the minor on the right 
hand side. That is, on the left we have a principal minor of B of size n + 1, and the principal minor 
of A on the right has size n — 1. Let X and Y denote the submatrices corresponding to these minors, 
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respectively, then we have det(X) = det(Y). Since Y is of Cauchy-type, the relation det(Y) = V c 
immediately follows. On the other hand, the computation of det(X) requires a longer preparation. 
To this end, we introduce the Cauchy-type (n + 1) x (n + 1) matrix $ with entries 

h a h b hgfc J c h b —|f|2 (A7 \ 

where a, 6 G N n . Recalling ([32]) . we see that 

X = $ + e(l + 2x c )- 1 e Cjn+ i + e(l - 2x c )~ 1 e n+ i iC , (48) 
i.e., X is a rank two perturbation of <£. Therefore, the determinant formula ()A.8p yields 

det(X) = det(d>) + e [j^2x~ c + TZT^ J + e (1 - 4xg) det(d>) ' (49) 

where the Cki's now denote the cofactors of <£. Using the special Cauchy-type form of we obtain 

det($) = -4a£(l - 4x2)- 1 P c | Wc z c | 2 , C c , n+1 = -(1 - 2x c )~ l V c Zo c z c , (50) 
together with the relations 



c c , c = p c i/ c i 2 n n[ Xc \ Xa) L Xc Xa) v ^ 5i ) 

„ = i (1 + x c + x a )(l - x c - x a ) 



/" T> II, |2 TT ( Xc Xa )( x c + ^a) , ko n 

It immediately follows that the determinant of X has the form 

det(X) = -(1 - Ax 2 c )~ 1 V c {Ax 2 c \lo c z c \ 2 + e{u c z c + ZJ c z c ) + e 2 ) . (53) 
Finally, by putting these formulae together, we end up with the quadratic equation 

Ax 2 c \oj c z c \ 2 + e(u c z c + uJ c z c ) + e 2 + 1 - Ax 2 . = 0. (54) 
Next, by solving equations (|45p and (|54p for z c , we find the following two solutions 



In order to select the right one, we proceed as follows. In the phase space region where the particles 
are far from each other, i.e., q\ S> • • • S> q n S> 0, moving with high relative momenta, i.e., p\ S> . . . S> 
p n 2> 0, the matrix L (|2ip is almost diagonal, therefore the diagonalizing matrix C7 (|26p is also nearly 
diagonal. Recalling equations ([30]) and (f34"p . we see that in the given phase space region the value of 
z c is very close to —1. Since z c is a smooth function over the connected phase space T*c, by invoking 
a standard continuity argument, we conclude that 

a=l 
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Having determined z c , we can find the form of the components f c and h c ()34p . too. As we discussed 
at the end of the previous subsection, by the non-uniqueness of U (|26|) , the non-zero complex quantities 
f c and h c are determined only up to a common phase factor ([36]). Therefore, purely for convenience, 
we may and shall assume that f c > for any c 6 N n . So, we can write 



i 



f c = e c \z c \? and h c = e c z c \z c \ 2 (57) 

with some 9 C E R. Combining this parametrization with (|32p . we obtain the following description of 
the matrix A. 

Lemma 2. With the aid of the X-dependent functions z c (56]). the matrix entries of A take the form 

A , _ Ja+6 b \ || 2l ff J _ -g a -0 h ^ 2l ff /rox 

2i# + A a -A 6 | ZaZb |a 2ig - A a + A b 

-tt-a,n+b — -tt-n+b,a — e Z[,\Z a Z h -f . O a fc, ^dyj 

215 + A a + A b 15 + A a 

where a, 6 E N n . 

Remark. Due to the presence of the 52-dependent second term on the right-hand side of equation 
(|59p. the matrix A can be seen as a C n -type non-trivial deformation of the usual Cauchy matrices. 
Spelling out the relation det(-A) = 1, the resulting Cauchy- type determinant formula might be of 
interest in other branches of mathematics and physics as well. 

3.3 Asymptotic phases and momenta 

Now let q(t) = (qi(t), . . . , q n (t)) E c be an arbitrary solution of the hyperbolic C n Sutherland dynamics, 
then by (|23p we can write 

{e 2 ^\ e^W, e-^W, . . . , e~ 2 ^} = a(Ae 2tL ), (60) 

where both L and A are computed at time t = 0. Since the diagonal entries of L = diag(A, —A) are 
not in decreasing order, we conjugate it by the unitary N x N matrix W := diag(l n , 7Z n ), where 7Z n 
is the unitary n x n matrix with entries (lZ n )a,b '■= $a+b,n+i ( a , b E N n ). Upon setting 

L := WLW~ l and A := WAW' 1 , (61) 

we see that the diagonal entries of L = diag(Ai, . . . , A n , — A n , . . . , — Ai) are decreasing, and for any 
a, b E N n we have 

Kb = Kb = fa{l + x a -x b )- l J b . (62) 

At this point Ruijsenaars' theorem [lj on the temporal asymptotics of exponential matrix flows is 
directly applicable. Recalling (|60p . for any c E N n we obtain the asymptotic relation 

e 2qc{t) ~ m c e 2tXc = e M^c)+2tx c ( t ^ QO ) ) ( 63 ) 

where the m c 's stand for the quotients of the consecutive leading principal minors of A. Since the cth 
leading principal minor of A has the form 

c 

V c :=ll\M 2 [J (l-(x a -x b )- 2 y\ (64) 

d=l l<a<b<c 
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we get 

c-l 



m 



_ V ^ _ \ f 12TT fi („ _ 

c 



a=l 



Therefore, for t — > oo we have g c (i) ~ gj" + tpj with asymptotic phases and momenta 

c-l 

2 — v " w 2 ~~ v "' 1 - 1 ' 2 



1 1 1 c_1 

g+ = -ln(m c ) = -ln(|/ c | 2 )--^ln(l-(x c -x a )- 2 ) and p+ = A c . (66) 



a=l 



By conjugating both L and A with the unitary N x N matrix TZn, we see that the diagonal entries 
of IZnLIZJ^ are in strictly increasing order, so the t — > —oo asymptotics can be handled similarly to 
the t oo case. It turns out that for t —¥ —oo we have q c (t) ~ q^ + tp~ , where 



1 1 c-l 

9c" = oMI^cl 2 ) - 2^2 ]n ( 1 ~ (. x c~ x a)~ 2 ) and p~ = -A c . (67) 

a=l 

Now we are in a position to formulate the main result of the paper. Indeed, recalling the form of 
z c ([55]) . the comparison of ([55]) and ([57]) immediately leads to the precise relationships between the 
asymptotic phases and momenta. 

Theorem 3. Take an arbitrary solution q{t) G c (t G K) of the hyperbolic C n Sutherland dynamics. 
For \t\ — )■ oo the particles move asymptotically freely, i.e., for any c G N n we have the asymptotics 

q c {t)~qt+tpt (t->±oo). (68) 
TTie asymptotic momenta satisfy the relations 

Pt = ~Pc and Pi > ■ ■ ■ > Pi > °> (69) 
and for the asymptotic phases we have 

c— 1 n n 

li = ~Qc ~^2 S (P7 -Pa>9)+ 6 (Pc ~Pa,9)+ d(p~ +p-,g) + 5(2p~,g 2 ) (70) 

a=l a=c+l a=l 

(a^c) 

with the 2-particle phase shift function 5(p,fi) = 2~ 1 ln(l + 4/i 2 p~ 2 ). 

4 Discussion 

In this paper we examined the scattering properties of the hyperbolic C n Sutherland model. Under 
certain technical conditions the classical repulsive particle systems are Liouville integrable, having only 
scattering states. Indeed, in general the asymptotic momenta provide sufficiently many independent 
first integrals in involution. However, the asymptotic momenta of the hyperbolic C n Sutherland model 
satisfy also the peculiar algebraic conditions p£ = —p~ (c G N n ). Following Ruijsenaars' terminology 
[5J, this distinguishing feature gives grounds for calling the C n -type model a pure soliton system. 
As in the ^4 n -type models, we expect that this stronger notion of integrability is responsible for the 
factorized form of the scattering map. 
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We find it interesting that the application of elementary linear algebraic techniques succeeds in 
revealing the scattering behavior of the C n model. Since the hyperbolic BC n Sutherland model with 
three independent coupling constants is also closely tied with the matrix Lie group U(n,n) (see [9]), 
we believe that the algebraic machinery presented in this paper can be extended to understand the 
scattering properties of the BC n model, too. Nevertheless, a complete classification of the pure soliton 
systems associated with the BC n root system appears to be a more challenging analytic problem. 
For further motivation we mention that these finite dimensional many-particle systems are closely 
connected with integrable field theories. For example, it is well-known that the family of A n -type 
Ruijsenaars-Schneider models describe the soliton solutions of the sine-Gordon equation (see e.g. |15j . 
|16j). Following [17J, it is also natural to speculate on the mathematically adequate description of the 
relationship between the scattering theory of the -BC n -type particle systems and the soliton dynamics 
of the boundary sine-Gordon models, at both the classical and the quantum level. We wish to come 
back to these issues in later publications. 

To explore the scattering properties of the C n model we adapted Ruijsenaars' approach [3] to the 
C n root system. Thus, it is quite natural that the matrix A appearing in Lemma 2 can be interpreted 
as the Lax matrix of the rational C n RSvD model. Recall that matrix A depends on 2n real parameters 
A = (Ai, . . . , A n ) and 9 = (9i, . . . , 9 n ) with Ai > . . . > A n > 0. Regarding the parameters A c and 9 C as 
canonically conjugated positions and momenta, the natural Hamiltonian associated with the matrix 
A takes the form 

H(\,9) = ±tr(A + A- 1 ) 

Now let us observe that, up to an irrelevant additive constant, this function can be identified with the 
rational limit of van Diejen's Hamiltonian (see equation (6) in [12J) with two independent coupling 
parameters. To prove the canonicity of A c and 9 C and to show the duality between the hyperbolic 
C n Sutherland and the rational C n RSvD models, one could imitate Ruijsenaars' original work [3] on 
the systems of type A n . However, for the ^4 n -type models it has been shown recently [18] that the 
symplectic reduction framework provides the most efficient way to establish the duality relation. For 
completeness, we now briefly outline the reduction picture underlying the duality of the C n models. 
The compact Lie group U(n, n)+ naturally acts on the manifold U(n, by conjugations. The lift of 
this action to the cotangent bundle T*[/(n,n)_ is Hamiltonian, admitting an equivariant momentum 
map J: T*U(n,n)- — ) u(n, n)* + = u(n, n)+. In the reduction picture our functional equation f)25[) 
corresponds to the momentum map constraint J = £ with the special choice £ f)15|) . Solving this 
constraint in a gauge in which A is diagonal, the outcome of the reduction procedure is the usual phase 
space of the Sutherland model (see e.g. [6], [7], [E]), and the Sutherland dynamics ([TJ is induced by the 
Hamiltonian H = tr(L 2 )/4. Note, however, that in Section 3 we solved the momentum map constraint 
in a gauge in which L is diagonal. As a result, the reduced phase space can be parametrized by the A c 's 
and the 9 c 's, and the dynamics of our interest (|7ip is induced by the Hamiltonian H = tr(A + A~ l ) /A. 
Therefore, by using two different gauges, one in which A is diagonal, and one in which L is diagonal, we 
obtain two different realizations of the reduced phase space T*U (n, n)_//,e U(n, n)+, whence a natural 
symplectomorphism between the phase spaces of the hyperbolic C n Sutherland and the rational C n 
RSvD models comes for free. The details, together with the generalization to the BC n systems, will 
be published elsewhere. 
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Notice that the action-angle duality between the Sutherland and the RSvD models has impor- 
tant consequences from the perspective of scattering theory, too. Based on the duality relation the 
scattering theory of the rational BC n RSvD models could also be understood. Moreover, as in the 
^4n-type models, the duality could greatly simplify the verification of the simplecticity of the M0ller 
wave transformations both for the Sutherland and the RSvD models, thereby the understanding of 
their scattering theory would be complete. 



A Appendix 

In this appendix we summarize some linear algebraic facts used throughout the paper. The proofs 
and further details can be found e.g. in [19J. 
For an n x n matrix X, let 

*(kk X) w 

denote the minor determinant of the pxp submatrix of X lying on the intersection of rows k% , &2 > • • • > 
k p with columns l\, I2, ■ ■ ., l p - Jacobi's theorem claims that there are simple relationships between the 
minors of X and the minors of its inverse, as described below. 

Theorem Al. Let X be an invertible n x n matrix, Y := and choose a permutation 

o={ kl k " "' kn )eS N (A.2) 



^ h h ••• I. . 

of the pairwise distinct indices k%, k^-, ■ ■ ■ , k n G N ra . Then for any p € {0, 1, . . . , n} we have 

Y ( h k 2 ■■■ k p \ _ sgn(q) x / k p+ i k p+2 ■■■ k n \ 
V h h ■ ■ ■ l P J det(X) \ l p +i l p +2 l n J ' 

where sgn(o") denotes the sign of permutation a. 

In this paper we frequently encounter Cauchy matrices and their perturbations. Recall that for the 
determinants of the Cauchy matrices we have 

1 \ n*<j(&-6)0»-*7*) 



where the ^'s and the 77^ 's are arbitrary complex numbers. The key formula that allows us to compute 
effectively the determinants of perturbed matrices is given in the following 

Theorem A2. Let X G C nxn be an invertible matrix and V, W G <C nxk be arbitrary matrices, then 
we have 

det(X + VW*) = det(X) det(l fc + W*X~ 1 V). (A.5) 



In particular, if X is perturbed by a multiple of the elementary matrix e a ^, then we can write 

det(X + ae a , fe ) = det(X) + aC a:b , (A.6) 
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where C a , b is the cofactor of X associated with entry X ajb , i.e., it is (— l) a+ times the (n — 1) x (n— 1) 
minor obtained by deleting the ath row and the 6th column of X. For analogue rank two perturbations 
we have 

det(X + ae a:b + 0e c , d ) = det(X) + aC a:b + (3C c4 + a(3(C a:b C c4 - C a , d C c , b ) det{X)-\ (A.7) 

Finally, if X is an invertible Hermitian matrix, the above formula simplifies to 

det(A + ae a , b + ae b>a ) = det(X) + aC a , b + aC a , b + |a| 2 (|C ib | 2 - C , a C b , b ) det(A)- 1 . (A.8) 

Since any quadratic submatrix of a Cauchy matrix is Cauchy again, the above determinant formulae 
offer a relatively painless way to compute the determinants of perturbed Cauchy matrices. 
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